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ABSTRACT 



We present an analysis of Spitzer IRS spectroscopy of 83 active galaxies from 
the extended 12/zm sample. We find rank correlations between several tracers of 
star formation which suggest that (1) the PAH feature is a reliable tracer of star 
formation, (2) there is a significant contribution to the heating of the cool dust 
by stars, (3) the H 2 emission is also primarily excited by star formation. 

The 55-90 vs. 20-30 spectral index plot is also a diagnostic of the relative 
contribution of Starburst to AGN. We see there is a large change in spectral 
index across the sample: Aa ~ 3 for both indices. Thus, the contribution to 
the IR spectrum from the AGN and starburst components can be comparable in 
magnitude but the relative contribution also varies widely across the sample. 

We find rank correlations between several AGN tracers. We find correlations 
of the ratios [Offl]A5007/[OIV]26/mi and [OIII]A5007/[NeV]14/im with the Sil 
strength which we adopt as an orientation indicator. This suggests that some 
of the [OIII]A5007 emiss ion in these Seyfe rts is subject to orientation dependent 
obscuration as found by lHaas et al.l (120051 ) for radio galaxies and quasars. There 
is no correlation of [NeV] EW with the Sil 10/im strength indicating that the 
[NeV] emission is not strongly orientation dependent. This suggests that the 
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obscuring materia l (e.g., torus) is not ver y optically thick at 14/iin consistent 



with the results of iBuchanan et al. I (120061 ) 



We search for correlations between AGN and Starburst tracers and we con- 
clude that the AGN and Starburst tracers are not correlated. This is consistent 
with our conclusion that the relative strength of the AGN and Starburst com- 
ponents varies widely across the sample. Thus, there is no simple link between 
AGN fueling and Black Hole Growth and star formation in these galaxies. 

The density diagnostic [NeV] 14/24 fim and [SIII] 18/33 jum line ratios are 
consistent with the gas being near the low density limit, i.e., ~ 10 3 cm" 3 for 
[NeV] and n e ~ few hundred cm -3 for [SIII]. 

The distribution of Sil 10 nm and 18 y mt strengths is consistent with the 



Ine distribution ol oil 1U nm an d io u , 
clumpy torus models of Sirocky et al. (j2008l ). 



We find a rank correlation between the [NeV] 14/im line and the 6.7 /im 
continuum which may be due to an extended component of hot dust . 

The Sy 2's with a Hidden Broad Line Region (HBLR) have a higher ratio of 
AGN to Starburst contribution to the SED than Sy 2's without an HBLR. This 
may contribute to the detection of the HBLR in polarized light. The Sy 2's with 
an HBLR are more similar to the Sy l's than they are to the Sy 2's without an 
HBLR. 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: spiral - infrared: 
galaxies 



1. Introduction 



The basic paradigm for nuclear activity in galaxies, an accretion disc around a super 
massive black hole (SMBH), is well-established. However, we have at best a sketchy under- 
standing of how SMBH are formed and what role nuclear activity plays in galaxy evolution. 
As we learn more about the formation and evolution of normal galaxies, and better under- 
stand the physics of active galactic nuclei (AGNs), it is becoming clear that these processes 
are intimately connected. The ubiquity of supe r massive black holes and the observed re- 
lation between black hole mass and bulge mass (IMerritt &: Ferraresdl200ll : iGebhardt et al. 
20001 ) suggest that nuclear activity may be closely connected to star formation in all galaxies. 
Accretion onto a black hole, which powers AGNs, requires available fuel in the nuclear region. 
Gas i nfall is also likely to be ass ociated with bursts of star formation (IHernquist fc Mihos 



19951 : iMihos fc Hernquist I Il996l ). Furthermore, AGNs deposit energy in the surrounding 



interstellar medium (ISM) via ionizing radiation, radio jets and perhaps winds, possibly 
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triggering star formation. Circumnuclear star formation has been observed in AGN and a 
link between nuclear activity and star formation is also supported by a recent finding that 
the mean stellar age of narrow- line AGNs is much younger than in similar normal spiral 
galaxies (iKauffmann et al.ll2003l ). However, the interaction between active nuclei and their 
host galaxies remains poorly understood. 

The mid-IR includes important diagnostics such as the PAHs , Sil features, and fine struc - 



ture lines (ISpinoriio & Malkan Ill992l ; IVoit Ill992bl ; iGenzel et al.lll998l : lLaurent et all 120001 1 . 



Previous mid-IR spectroscopic studies of AGN have shown that these diagnostics allow sep- 



obscuration of the central torus (Genzel et al. 1998; 


^aurent et al. 2000; Clavel et al. 


2000 


Sturm et al. 2002; Haas et a 


2005; 


Weedman et al. 


2005: Ode et al. 


20061: Shi et al. 


2006 


Hao et al. 2007; Spoon et al. 


2007). 



This paper is part of a series on a Spitzer study of the Extended 12/zm Sample of Seyfert 
galaxies. iBuchanan et al. I (120061 ) presented low resolution IRS spectra of the first 51 (of 83) 
sources and demonstrated the variety of mid-IR spectral shapes due to the relative mix of 
comparable contributions from starburst and AGN. They examined the ratio of IR to nuclear 
radio flux density and showed that the ratio is a factor of two larger in Sy l's than Sy 2's 
between 15 and 35 /xm suggesting the torus emission is anisotropic by a factor of two at 



these wavelengths. Additional preliminary results have been presented by IBuchanan et al 
(120081). Other papers curren tly planned in the series will present (1) an atlas of all the 



data (IGallimore et al. Il2009al ). (2) models of the mid- to far - IR SEDs which combine both 
starburst and AGN contributions (IGallimore et al. I l2009bl ). and (3) UKIRT observations 
of the near-IR emission and discuss constraints on the relative mix of starburst and AGN 



contributions to the spectrum (IBuchanan et al. 1120091 ). In this paper we examine diagnostic 
features (PAH, Sil, fine structure lines, H 2 , continuum slope) in the IRS spectra and discuss 
the implications for the AGN-Starburst connection and obscuration of the cen tral engine. We 
note t hat independent groups have also pub 
(120071 ) used observations of some Sy l's from 



ished results based on our data. Maiolino et al. 



Buchanan et al. I (120061 ) as a comparison sample 



for their study of high luminosity QSOs. IDeo et al.l (120071 ) obtained low resolution Spitzer 
IRS observations of 1 2 Sy 1.8 and 1.9 and c ombi ned the results w i th ar chival data of 57 
Sy l's and Sy 2's from IBuchanan et al. I (120061 ) and lWeedman et al. 1(l2005h. High resolutio n 



IRS obse rvations of a sub sample of 29 sources are presented by iTommasin et al. I (120081 ). 
Recently IWu et al. I (1200 9) published a study of 103 Seyferts which included our IRS data. 



Gallimore et al. I (j2009al ) show that there are systematic d ifferences between our values for 



the PAH and Sil features and those of IWu et al. I (120091 ) which are due to the different 



analysis techniques used (PAHFIT vs. Spline-fitting). Nevertheless, we agree on the overall 
trends which are common to both papers. 
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2. Spitzer Observations 



2.1. The Sample of AGN 



The properties of the sample are discussed by lBuchanan et al. I (120061 ) and lGallimore et al 



(j2009al ). The sample co mprises all g a laxies identified as Seyferts or LINERS from the ex- 
tended 12 /im sample of iRush et al.l (119931 ) that have cz < 10000 km s _1 . The potential 
advantage o f the 12 um sample is its perceived resistance to wavelength-dependent selec- 
tion effects (jSpinoglio fc Malkanlll989l ). Note that there may in fa ct be selection effects i n 



the sample due to anisotropic emission from the torus at 12 /im ( iBuchanan et al. II2006I ). 



The redshift criterion was included to ensure that the smallest possible physical size was 
probed by the fixed Spitzer apertures in order to better exclude host galaxy emission. The 
data include IRAC imag ing, IRS spectral mapping, and MIPS SED spectra as described by 



Gallimore et al. I (I2009al ). 



The spectral class ifications were ob t ained b y a thorough search of the literature for each 
galaxy; see Table 1 of iGallimore et al. I (l2009ah . For the purposes of comparing the Seyfert 
types we consider Seyfert types l.n (i.e. types 1.0, 1.2, 1.5, 1.8, and 1.9 as well as narrow 
line Seyfert l's) to be Sy l's, and objects with no broad permitted lines to be Sy 2's. We 
consider the properties of the Seyferts with hidden broad line regions (HBLRs) separately 
and compare them to the Sy l's and Sy 2's. We refer to the sum of the Sy l's, Sy 2's, and 
hidden Sy l's as the combined sample. The resulting sample contains 83 sources of which 7 
are now classified as HII galaxies, 6 are LINERS and the remaining 70 are Seyferts (27 Sy 
1-1.5, 17 S y 2 with HBLR, 20 Sy 2 without HBLR, and 6 Sy 1.8-1.9. The sample is listed in 
Table 1 of IGallimore et al. I (l2009ah . 



2.2. IRS Observations and Data Reduction 



The sample galaxies were observed by Spitzer IRS (IHouck et al.l 120041 ) in spectral map- 
ping mode, using the Short-Low and Long-Low modules. The modules cover a wavelength 
range of ~5 - 36 /im, with a resolving power (A/ A A) ranging from 64 to 128. The integration 
times were 6 seconds pe r slit pointing. Details of the observing and reduction are given by 
Gallimore et al. I (j2009al ) and are summarized here. 



The spectral maps were constructed to span > 10" in the cross-slit direction and centered 
on the target source. The observations were stepped perpendicular to the slit by half slit- 
width spacings. For the SL data, the mapping involved 13 observations stepped by 1.8" 
perpendicular to the slit, and for the LL data, 5 observations stepped by 5.25" perpendicular 
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to the slit. The resulting spectral cubes span roughly 25.3" x 54.6" (5.3 - 14.2 fim) for the 
SL data and 29.1" x 151" (14.2 - 36 /jm) for the LL data. 

The raw data were processed through the Spitzer BCG pipeline, version S15.3.0. Sky 
subtraction employed off-source orders. In SL and LL observations, the source is centered in 
first or second order separately, with the "off-order" observing the sky at an offset position 
along the slit orientation. Sky frames were constructed using median combinations of the 
off-source data and subtracted from the on-source data of matching order. No detectable 
contamination appeared in the sky frame images. 

Data cubes were constructed and edited for cosmic-rays and spurious pixels using the 



CUBISM software package developed by the SINGS IRS team (ISmith et al.ll2007bl ). We used 
synthetic, 20" diameter circular apertures to extract the spectra. 

To calibrate the flux scale, we employed archival IRS staring and spectral mapping ob- 
servations of the calibrator stars HR7341 (SL) and HR6606 (LL). The staring mode spectra 
were extracted using SPICE (Spitzer Science Center 2006), and the spectral mapping obser- 
vations were processed as described above. The flux calibration curve was determined by the 
ratio of the flux-calibrated staring mode spectra to the uncalibrated cube-extracted spectra. 
This flux calibration was applied to all target sources observed in spectral mapping mode. 
Based on inspection of the spectral overlaps between modules and orders the calibration is 
self-consistent to within 5% across the IRS spectrum. The only exception was NGC 4151, 
SL order 1, which had to be scaled down by ~ 15% to match the neighboring SL order 2 
and LL order 2 spectra, perhaps owing to a residual pipeline error. After applying a uniform 
scaling to SL order 1, the spectrum matched well the neighboring spectral overlaps. 

Archival data of our sample, such as obtained by the SINGS collaboration, were pro- 
cessed in the same manner. There are however a few sources for which only staring mode 
observations are available for the SL module: NGC 526A, NGC 3227, NGC 4941, IC 5063, 
NGC 7172, and NGC 7314. The staring mode spectra were extracted using S PICE and scaled 



to mat ch IRAC 5.8 and 8.0 /zm photometry obtained as part of this project (iGallimore et al. 



2009ah 



2.3. MIPS-SED Observations and Data Reduction 

The sample galaxies were also observed using the MIPS camera in SED mode, which 
provides very low resolution (R ~ 20) slit spectroscopy covering ~ A55-90 /im. The slit- 
width for SED mode observations is 20". Integration times ranged from 3 to 10 seconds 
depending on the IRAS 60/um flux density. Observations include six pairs of on-source and 
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off-source positions (1-3' away from source) for background subtraction. 



Basic data processing, including background subtraction and flux and wavelength cali- 
bration, was handled by the Spitzer pipeline (version S.14.4). The final, low resolution spec- 
tra were extracted from the long-slit images using a three-pixel (~ 30") synthetic aperture 
perpendicular to the slit. For each source, the position of the galaxy nucleus was identified 
first manually and then by a Gaussian fit to determine the center of the extraction aperture 
to a fraction of a pixel. 

Note that flux calibration is derived from observations of standard stars and so is most 
accurate for compact far-infrared sources. Formally, the accuracy for point sources is better 
than 10 %, but degrades to ~ 15% for extended sources with galaxy-like surface brightness 
profiles JLu et al. Ibooi l. 



2.4. PAHFIT Spectral Decomposition 



We used the PAHFIT spectral decomposition code (ISmith et al.ll2007al ) to fit and ex- 
tract line features of the Seyfert IRS spectra. The publicly available release of PAHFIT 
permits silicate features only in absorption. Silicate 10 /im and 18 /im emission commonly 
appears among the sample, however, and we therefore modified PAHFIT to include a warm, 
optically thin (r(10 /im < 1) dusty medium. To accommodate the range of observed Sil 
emission profiles, the warm dust model includes two components: (1) an optically thin warm 
medium, characterized by a single temperature and 10 /im opacity, and (2) a cold, partially 
covering dusty medium parameterized solely by covering fraction (i.e., the cold dust is mod- 
eled to produce no significant emission over the IR S band). Bo t h the c old and warm dust 



components obey the extinction law described by ISmith et al.l (l2007al ). Our goal was to 
reproduce the observed silicate emission profile shapes without regard to determining the 
physics of the warm dust (other than its presence and emission strength), and so, for simplic- 
ity, the spectrum of the warm dust, modified by partial foreground absorber, was calculated 
assuming a simple slab geometry. 

Except for this added warm dust component, we used PAHFIT with its default settings. 
PAHFIT robustly decomposes the spectrum into multiple temperature dust components, 
silicate features, emission lines, and PAH dust features. The decomposition proceeds using a 
standard non-linear fitting algorithm. PAHFIT was specifically designed to analyze moderate 
resolution IRS spectra and carefully limits parameters, such as line widths and centroids, 
that might drift to unphysical values during the fit. 



The final data products include tables of line strengths and equivalent widths. We 
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further modified PAHFIT to produce line-subtracted continua. These continua were used to 
calculate silicate strengths and spectral indices. In addition, we also modified the PAHFIT 
output to produce 3cr uppe r limits for weak or undetected lines. These data tables are 
provided in iGallimore et al. I (j2009al ) . Note that the line and feature fluxes are not corrected 
for extinction. 

Table [1] summarizes the results of the PAHFIT decomposition, giving the median values 
of a range of properties for the continuum emission, Sil feature, star formation tracers, AGN 
tracers, and diagnostic lin e ratios. The median value s were calculated based on the compu- 
tation of a survival curve (IFeigelson fc Nelson Ill985l ). As the line ratios include both upper 
and lo wer limits , we e mployed the formalism for interval-censoring given by iFinkelstein 
( 119861 ) and ISun I (119961 ). as implemented in the interval package by M. P. Fay & P. A. Shaw 



(http: / / cran.r-project.org/web/packages/interval/index.html 



Tables EH9] give the results of statistical tests for the significance of differences in the 
distribution of those properties for the different Seyfert subgroups. These tests comprise 
Kolmogorov-Smirnov two-sample tests for uncensored data, or log-rank tests for interval- 
censored data. 

Table [10] gives the significance of rank correlations between the properties for the com- 
bined sample. We used Kendall's r ank correlation statistic adapted for interval-censored data 
( llsobe. Feigelson. fc Nelson I Il986l ). reported as a z-score and two-sided significance level. 



In this paper, we search extensively for patterns in the data by exploring bivariate cor- 
relations among measurements: for example, PAH eqw vs. 20-30 /zm spectral index. There 
are well-known parametr ic and non-parametric statistical measures of correlation (see, e.g., 
Sprent fc Smeeton I l2000h . Parametric methods specifically apply to cases where the data 
are believed to be related through some known and prescribed functional form; for example, 
the well-known Pearson correlation coefficient measures the underlying linear relationship 
between measurements. One however expects a fairly complex relationship, if at all, between 
a spectral index and an eqw, and so such parametric statistics would not be expected to ap- 
ply in the present analyses. For example, data that are tightly but n on-linearly correlated 
will likely show low linear correlation measurements (lAnscombe Ill973l ). 



In contrast, rank correlation coefficients provide a robust, non-parametric measure of 
association between measurements. They are better suited for the present study, as they 
measure simply the degree to which one measurement increases (or decreases) as another 
increases, independent of the degree or functional form of the underl ying association . To 
our knowledge, only Kendall's rank correlation statistic (Kendall's r; ( Kendall 19381 )) has 
been generalized to account both for upper limits (non-detections) and lower limits (e.g., line 
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ratios where the divisor is an upper limit) in one or both measurements and so is ideal for this 
study. We have i mplemented the generalized Kend all's rank correlation statistic following 
the formalism of llsobe. Feigelson. Sz Nelson I Jl986h . In the results tables, we report the 
standard normal score (z) of the statistic as well as the probability (p) that score would have 
occurred by chance (i.e., the probability of no correlation). We chose to be conservative in 
the interpretation of the correlation statistics and adopt the threshold p < 1% to indicate 
evidence of rank correlation. 



3. RESULTS 



3.1. The Starburst AGN Connection 



3.1.1. PAH Emission as a diagnostic of Star Formation 

The PAH emission is thought to be produced when small grains (containing ~ 50 Carbon 
atoms) are heated by singl e UV photons from hot stars and radiate though their vibrational 
modes ( Tielens et al. 1999). It is cl ear that PAH emissio n is a tracer of star forming regions 



in our Galaxy ( Tie 



2004; Brandl et al. 



ens 



et all Il999h and other galaxies (jLaurent et al.l 120001; iPeeters et al.l 



2006 



Dale et al.ll2006l : iDraine et al.ll2007l : ISmith et al.ll2007ah The PAH 



luminosity is proportional to the star formation rate in star forming galaxies (IRoussel et al. 



2001 




3randl et al. 


2006; 


(1992a 


) argues that PAH 



Calzetti et al.l 120071 : IShi et all 120071 : iFarrah et al.l 120071 ). 



Voit 



destroyed, and the PAHs in AGN can only survive when they ar e shielded by high co lumns 



of X-ray absorbing gas, e.g., in the plane of the obscuring torus. ISmith et al.l (j2007al ) notes 
that the ratio of 7.7/11.3 micron PAH seems to be lower in some AGN than in starburst 
galaxies. Since the 11.3 micron PAH are produced by larger grains than the 7.7 micron PAH 



( IDraine fc Li 1120071 1 this suggests that the smaller grains are indeed being destroyed by the 
AGN. Thus, PAH emission seems to be a good indicator of star formation in galaxies with the 
caveat that it may be significantly suppressed in AGN. Siebenmorgen, Kriigel & Spoon (2004) 
compare radiative transfer models which include dust destruction for starbursts and AGN. 
In their models, some PAH emission is produced within clouds illuminated by an AGN, but 
the strength of the PAH is lower by a factor of about 5 than that produced in the starburst 
models with the same luminosity and A„. Finally, the equivalent width of PAH features 
has been shown to be useful in diagnostic diagram s as a discriminant between starbursts 



and AGN as the energy source for infrared emission ( Genzel et al.l Il998l ; lLaurent et al.l 12000 



Peeters et al.l I2004J : lArmus et al. 1 120071 ; ISpoon et al. 



20071 ). 



The 6.2 micron PAH feature is thought to be relatively "unconfused" by other spectral 
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features (e.g., Peeters et al. 2004) and so should be a good measure of the amount of star 
formation. We show the distribution of 6.2 /im PAH EW for the Sy l's, Sy 2's and hidden 
Sy l's in Figure [TJ We see that the Sy 2s show a distribution to higher values of PAH 
EW. We find that the distributi on of PAH EW i n the Sy l's and Sy 2's are significantly 
different (Table |2J) in contrast to IWu et al. I (120091 ) who find no significant difference in the 
E Ws. This might be due to the systematic differences in the measured PAH EWs discussed 
by iGallimore et al. I (j2009al ). The Sy 2's which show the high EW are those without an 
HBLR; the Sy 2's with an HBLR have a distribution which is similar to that of the Sy l's 
(Tables CD, [2|3|) . The higher values of PAH EW for the Sy 2s seems mainly due to the lower 
continuum fluxes in the Sy 2's relative to the Sy l's since the PAH luminosities (Figure [2]) 
are similar in the Sy 2's and the Sy l's. This result is consistent with the findings of Clavel 
et al. (2000) for a sample of Seyferts studied using ISO. Figures [J and El also show that the 
Seyferts show a similar distribution in 6.2 /im PAH luminosity as the iBrandl et al. I (120061 ) 
starbursts. However, the Seyferts show lower 6.2 /im PAH EW due to the higher continuum 



contribution from the AGN. These results are consistent with the results of Buchanan et al. 



( 120061 ) that principle component analysis required a higher star formation contribution to 
the SED of Sy 2's. 



3.1.2. The PAH 6.2/11.3 ratio 



Th e 6.2 jtrni PAH featu re is produced by grains with smaller sizes than the 11.3 /im PAH 
feature (jDraine fc Li 1 120071 ). The smaller grains are more quickly destroyed by exposure to 
the AGN spectrum. T hus, the ratio of 6.2/ 1 1.3 PAH should b e an indicator of the destruction 
of grain s by the AGN dSmith et al.ll2007af) . IWu et al. I (120091 ) show that their Seyfert sample 
and the lBrandl et al. I (120061 ) sample of starburst galaxies have similar average values for the 
11.3/6.2 PAH ratio. 

We calculate the 6.2/11.3 PAH ratio using the fluxes of the features. The distribution 
of the 6.2/11.3 PAH ratio is consistent with being the same in the different subgroups. In 
Figure [3] we plot the ratio of the 6.2 and 11.3 /im PAH strength vs. the [NeV] luminosity 
which traces the strength of the AGN ionizing continuum (§ I3.1.6p . We see no correlation of 
the PAH ratio with the [NeV] luminosity. In FigureHJwe plot the 6.2/11.3 PAH ratio against 
the Sil 10 /im Strength (see §3.2.ip . We do not see a correlation between the PAH ratio and 
the Sil 10 /im Strength (Table [TO]) . This suggests that the PAH emission that we see in these 
Seyferts is not orientation dependent. Any effect of the AGN on the distribution of PAH 
grains must be confined to a small volume so that the integrated measurements studied here 
are not affected. 
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3.1.3. PAH Associated with Emission from Cool Dust 



Schweitzer et al.l (120061 ) show that QSOs and ULIRGs lie on the same correlation be- 
tween 7.7 /im PAH flux and 60 /im continuum flux, suggesting that the FIR is dominated by 
dust heated by star formation. In Figure [5] we show a strong rank correlation between 6.2 
/im PAH luminosity and the integrated luminosity (from IRAS) at 60 /im for our sample. 
We also find a rank correlation between the flux of the 6.2 /im PAH feature and the 60 /im 
continuum (Table [TO]) . 

We also show the 6.2 /im PAH EW vs. the spectral index of the MIPS SED 55-90 /im 
in Figure El The slope of the MIPS SED indicates the relative importance of cold vs. warm 
dust to the FIR SED. Positive slope is a red MFIR SED. The starburst galaxies all lie at 
large PAH EW. In the combined sample, we see a rank correlation between PAH EW and 
MIPS slope, with higher PAH strength (more star formation) proportional to more positive 
(redder) MIPS slope (more emission from cool dust). These two correlations (Figures and 
E]) are consistent with a significant contribution from stars to the heating of the cool dust 
which produces the FIR emission in our sample. 



3.1.4. PAH EW vs. 20-30 \im spectral index: a measure of the relative strength of the 

starburst and AGN components 

20 microns is just r edward of the "knee bend" we see in the SED of many Sy Is and 2s 
( IBuchanan et al. II2006I ). The 20-30 micron spectral index a 20 _3o, was determined using a 
power law fit to the line-subtracted continuum between 20 and 30 microns rest frame. This 
spectral index characterizes the relative contribution of warm a nd cool dust to the mid-IR. 
Note that the values of a 20 _3o used in IBuchanan et al. I (120061 ) are actually a 2 o-30 + 2.0. 
The correct values of a 2 o-30 are used here. The s imilar 15-30 (im spectral index has als o 
been found to be an effecti ve AGN discrimin ator (IBrandl et al. 1 120061 ; ICleary et al.l 120071 ). 
Using a subset of our data, iDeo et al.l (120 07') show that there is a correlation between the 
6.2 /im PAH EW and a 2 o-3o- IWu et al. I (120091 ) also find a correlation between the the 11.3 
/im PAH EW and the 15-30/im spectral index. We plot the 6.2 /im PAH EW vs. a 2 o_3o in 
Figure We see the different subsamples tend to lie in different parts of the diagram. The 
Sy l's and Sy 2's with an HBLR generally have a 2 o-3o < —0.5 and cover a large range in 6.2 
/im PAH EW. On the other hand the Sy 2's and Liners tend to have high values of 6.2 /im 
PAH EW and cover a large range in a 2 o-3o with the Sy 2's showing a strong rank correlation 
(p < 0.1%) between the two properties. We also see a rank correlation between 6.2 /im 
PAH Luminosity and a 2 o-3o (Figure [HJ). Thus we see rank correlations between the PAH 
EW and a 2 o-3o which are both measures of the relative strength of the starburst and AGN 
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components, and also PAH luminosity (a measure of the absolute strength of the starburst 
component) and «20-30- 

In Figure [7] the sources in the lower left are those with a low starburst to AGN ratio. 
Thus, the value of «2o-30 between -1 and -2 is likely t o repre sent the spectrum of the "torus", 
consistent with values calculated by iNenkova et al. I (120081 ) . 

In Figure [9] we plot the 55-90 /im spectral index vs. 20-30 /im spectral index. We see 
a correlation in the sense that objects which are redder in 20-30 tend to be redder in 55-90. 
This is consistent with the hypothesis that star formation contributes to both the 20-30 and 
55-90 spectral regions. Thus, the 55-90 vs. 20-30 spectral index plot is also a diagnostic of 
the relative contribution of Starburst to AGN. We see there is a large range in spectral index 
across the sample: Aa ~ 3 for both indexes. Thus, the contribution to the IR spectrum 
from the AGN and starburst components can be compar able in magnitude but the relative 



contr i bution also varie s widely across the sample (see also iBuchanan et al. 1 12006c iDeo et al. 



20071 : IWu et al. 1 120091 ) . 



3.1.5. The Relationship between Molecular Hydrogen and Star Formation 

Figure dD] shows the 6.2 /im PAH EW vs. the H 2 S(l) 17 /im EW. The H 2 S(l) line has 
higher EW in the Starbursts and Liners than the Seyferts. The distribution of the fluxes and 
luminosities of the ^5(1) 17 /im line are consistent with being the same in the subclasses 
of Seyferts (Tables | 



We see that the H 2 5(1) EW is proportional to the PAH EW. iTommasin et al. I (120081 1 
show a correlation between the luminosities of the 11.3 /im PAH and H 2 S(l) 17 /im lines in 
a subsample of these sources. We also find a rank correlation between the luminosities of the 
6.2 /im PAH and th e H? S(l) 17 /an l i ne (a rank correlation is also found in the fluxes; Table 



ITUT) . We note that ISchweitzer et al.l (120061 ) have also reported a correlation between PAH 
and H 2 in the higher luminosity QUEST QSOs. Since the PAH is a tracer of star formation 
this suggests that the H 2 line is also associated with star formation. 

In Figure[TT]we plot a 20 -3o vs. the H 2 S(l) 17 /im EW. We see a correlation in the sense 
that objects with redder continuum slopes have larger H 2 S(l) EW. Since red continua are 
associated with cold dust from stars, this correlation is also consistent with the molecular 
hydrogen line being associated with star formation. 
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3.1.6. AGN Tracers 



Following IVoit I Jl992bh and ISpinoglio fc Malkan I (Il992l ). we adopt the high ionization 
lines [NeV] (97.1 eV) and [OIV] (54.9 eV) as tracers of the AGN in these objects. The 
vast majority of stars do not produce photons suf ficiently energetic to excite [NeV] and 
[OIVl , with the ex c eptio n of young Wolf-Rayet stars ( Schaerer fc Stasihska 19991 ). However, 
Abel fc Satyapal I (120081 ) show that even in AGN with a strong starburst component, the 
observed ionization parameters constrain the [NeV] emission to be dominated by the AGN. 

It is also possible for [NeV] and [OIV] to be excited by shocks with velocities > 200 



km s (lAlen et al. 1 120081 ). However, as Allen et al. note, the predicted line ratios are 
not in agreement with observations. Detailed studies of th e IR spectrum of Seyferts are 



inste ad consistent with AGN photoionization of the NLR (IGroves. Dopita fc Sutherland 
20061 ). In addition, the line ratios [NeV] /[Nell] and [OIV] /[Nell] have proved empirically 



to provide good separation between AGN and Starbursts (IGenzel et al. 



1998; Sturm et al. 



20021 . 120061 ) . And the [OIV]26/im line has also shown to be correlated with other AGN 



tracers - the [OIH]A5007 line, radio emission, and t he 2-10 keV X-ray flux in other samples 
rtMelendez et alJbood : iDiamond-Stanic et al. Il2009h . 



In Figure [T2] we show a rank correlation between the luminosity of the [NeV] 14/zm line 
and the 8.4 GHz radio continuum. A rank correlation is also present between the fluxes 
of the [NeV] 14/im line and the 8.4 G Hz radio continuum (Table [TOl) . The 8.4 GHz radio 
emission is not orientation dependent (IThean et al. Il200ll ) and is an unbiased measure of 
the isotropic AGN luminosity. The rank correlation of the [NeV] 14 fim line with the 8 
GHz radio emission suggests that the [NeV] line emission is also a useful measure of AGN 
luminosity (see also § 13.2.3ft . 



In models of the obscuring torus, th e flux in the near-m i d IR (i.e., 



is due to hot dust heated by the AGN ( Pier fc Krolik 



1992 



Efstathiou fc Rowan- Robinson 111995; van Bemmel fc Dullemond 2003; Ncnkova et al. 



1 0/im) 



Granato fc Danese 



994 



2002 



20081 ). Thus, we adopt the 6.7 //m continuum flux as an AGN tracer. In § 13.2.41 we show a 
rank correlation between the 6.7 fim continuum and the [NeV] 14 /im line emission. 

In Figure [13] we see rank correlations between the luminosities of the two IR AGN trac- 



ers the [NeV] Uu m and [OIV] 26^m lines with the op t ical [OIII1A5007 



Gu et al. 



2006 



Kewlev et al. 



2 



ments taken from (lAxon 2009lbrivate communication!; iDeo et al 



2001 



Kim et al. 



1995l;lKirhakos fc Stei 



200 



7: 



emer 



ine ([QUI] measure- 



de Griip et al.lll992l: 



199 



0; 



Melendez et al 



20021 : iTranl I2OO3I : IVeron-Cettv fc Veron 



2008; 



Panessa fc Bassani 



Winkler 



1992 



Whittle 



2002 



1992 



Smith ct a 



Young et al. 



2001 



19961 )). The rank correlations are also present 



but weaker between the fluxes of the lines. This confirms that the [NeV] 14/im and [OIV] 



13 



26/im are useful AGN tracers. 



3.1.7. Comparison of Starburst and AGN Tracers 



Correlations between AGN aud starburst tracers have been fou nd for QSOs flHaas et al. 
20031 : ISchweitzer et al.lbood : iNetzer et al1l2007l : llVlaiolino et~aDl2007h and Seyferts JlVlaiolino eTal 



20071 ). Figure [T41 shows a rank correlation between the luminosities of the 6.2 /im PAH and 
the [NeV] 14 /im line in the combined sample. We also find rank correlations between the 
luminosities of (1) the 6.2 /im PAH and the [OIV] 26 /im line (2) the [NeV] 14 /im line 
and the 60 /im continuum, and (3) the 6.2 /tm PAH feature and the 6.7/mi continuum (not 
plotted, see Table [10]). These correlations seem to indicate that the starburst and AGN 
luminosities are correlated in these Seyferts. However, our statistical tests (Table [TO]) show 
that there are no significant correlations in the flux of these AGN and Starburst tracers. 
Figure (TS] shows that the fluxes of the 6.2 /im PAH and [NeV] 14 /im line scatter over an 
order of magnitude and are uncorrelated. So although AGN tracers are internally correlated 
in both flux and luminosity, and Starburst tracers are are internally correlated in both flux 
and luminosity, the fluxes of the AGN and Starburst tracers are not correlated. We con- 
clude that the correlations in Starburst and AGN luminosity are due to the dependence on 
distance and are not real. This is consistent with our results that the relative strength of the 
AGN and Starburst components in these Seyferts vary widely across the sample (§ I3.1.4p . 
The lack of a correlation between the AGN and star burst suggests that in these AGN there 
is no direct link between AGN fueling, black hole growth, and star formation. 



3.1.8. The [NeV] /[Nell] ratio vs. PAH EW Diagram 



Genzel et al.l (119981 ) showed that a diagram which plots the ratio of a high to low ioniza- 



tion li ne vs. a PAH EW is able to separate Starbursts from AGN from ULIRGs. iGenzel et al. 
Jl998h used [OIV]/NeII] vs. the 7.7/mi PAH EW. In Figure [U we plot 6.2/iin PAH EW vs. 



[NeV] 14/im/ [Nell] 12/im. We see that there are significant differences across our sample with 
objects with higher AGN to starburst ratio in the lower right and the objects with lower AGN 
to starburst ratio in the upper left. Note that the Liners have low [NeV]14/im/[NeII]12/im 
and tend to fall to the left of the other objects. 
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3.2. Orientation, Obscuration, and Constraints on the Torus 



3.2.1. 6.2 fxm PAH vs. Sil 10 pm Strength 
The strength of the Sil 10 am feature is potentially a powerful probe of the absorbing 



mate rial (IPier & Krolik I Il992l ; lLaor Draine I Il993l ; iLevenson et a l. 2007; llmanishi et al. 



20071 ). We adopt the definition of Sil strength given by ISpoon et all (120071 ). 



Si 



Sil 



In 



/obs(10) 



(1) 



/cent (10) 

where / Q b s is the observed flux in the feature and / con t is the interpolated continuum flux 
at 10 /im. Recent Spitzer res ults suggest that different types of objects show systematically 
different Sil 10 /im strength (jShi et al.ll2006l ; lHao et al.ll2007l ) with QSOs showing strong Sil 
emission, Sy Is showing mostly weak Sil emission, Sy 2s sh owing mostly weak to moderate 
Sil absorption, and ULIRGS showing strong Sil absorption. ISpoon et all (120071 ) have shown 
that objects fall mainly in two branches on the 6.2 /im PAH vs. Sil 10 /im Strength diagram. 

Figure [T7| shows the distribution o n the 6.2 um PAH vs. Sil 10 /im Strength diagram. 
Our results are consistent with those of ISpoon et all (120071 ). As expected, the Sy l's tend to 
show weak Sil emission, and the Sy 2's tend to show Sil absorption. The KS test (Table [2]) 
confirms that the distribution of Sil strength i s very diff e rent i n the Sy l's and Sy 2's. 



These results are in qua litative agreement with IWu et al. I (120091 ). though as discussed by 



Gallimore et al. I (j2009al ) there are systematic differences in the measured Sil strengths. The 
only significant difference between the Sy l's and the Sy 2's with an HBLR is the Sil strength 
(Table E]). 

Two of the HBLR Sy 2's show Sil emission - F01475-0740 and M-3-58-7. In these 
objects, the obscuring material (torus) may be oriented at an intermediate inclination. In 
addition, four objects classified as SI. 8-1. 9 show Sil absorption - NGC1194 (S1.9), M-2-40-4 
(Syl.9), NGC1365 (S1.8), and NGC5033 (SI. 8). These are also likely to be fairly edge-on 
(§ 13.31) . About half of the objects identified as Liners and Starbursts show weak Sil features, 



with strong absorption seen in MRK938 (HII), NGC3079 (LINER), NGC5005 (LINER), and 
strong emission seen in NGC4579 (LINER), NGC5494 (LINER), and CGCG381-051 (HII). 



3.2.2. Sil 10 um Strength vs. Sil 18 pm Strength 



Sirocky et al.1 (120081 ) showed that the observed relationship between the two Sil features 
18 /im and 10 /im can be a powerful probe of the properties of the obscuring material. Fig- 
ure [18] shows the distribution of the Sil strengths for the different subsamples of galaxies. 



- 15 - 



We see the Sy l's mainly to the upper right with emission in both Sil features, and the 
Sy 2's mainly to the lower left with absorption in both Sil features. In Figure [19] we com- 
pare our results for the two Sil Strength features with simple spherical obscuration models 
( ISirocky et al.ll2008l ). The spherical models have the advantage that they reduce the number 
of free parameters and still g ive good agreement with the more detailed clumpy torus models 
( jNenkova et al. 1120021 . 120081 ) . We see that the data lie below the smooth torus models and 
within the parameter space of the clumpy torus models. Our data are inconsistent with the 
smooth torus models but are consistent with the clumpy torus models. Within the context 
of the spherical clumpy models, there is a preference for a small number of clumps along the 
line of sight (N ~ 1 - 3). 



3.2.3. Density and Obscuration Diagnostics 



The [NeVl 14/24 and [SIII] 18/33 line ratios are density diagnostics (IVoit I Il992b 



Spinoglio & Malkan 


199 


2j), and 


tinction indicators P 


/oft 


1992bl: 



and in principle, in the low density limit ca n be used as ex 



will probe high excitation gas, most likely in the narrow line region, while the [SIII] can be 
excited by stars and may be produced on a larger spatial scale. Sturm et al (2002) examined 
the [NeV] 14/24 ratio in a sample of Seyferts using ISO data and also concluded that the 
[NeVl emitting gas in both Seyfert types was in the low density limit. On the other hand 
Dudik et al.l (120071 ) find that type 2 AGNs tend to have lower ratios than type 1 AGN, and 
some of the type 2 AGNs are below the low density limit possibly due to differential absorp- 
tion. Tommasin et al. (2008) using high resolution IRS spectra on a subsample of the 12 
jum sample, report [NeV] line ratios slightly above the low density limit, giving densities of 



n,, 



10 3 cm 



Our data (Figure 120]) and Tables CD - EJ) are consistent with Sy l's, Sy 2's and hidden 
Sy l's having the same distribution of these diagnostic ratios. Both line ratios seem to 
scatter around their low density values giving for [NeV] n e ~ 10 3 cm -3 and for [SIII], n e ~ 
few hundred cm -3 . We also see no evidence for a correlation between the values of these 
ratios (Table [TO]) . 



In Figure |2~T1 we plot the [NeV] 14/24 line ratio against the Sil 10 /im feature which is 
an indicator of absorption along the line of sight ( §3.2.11) . We see that there is no correlation 
between these two quantities (Table [TO]) . This suggests that extinction does not have a 
significant effect on the observed [NeV] 14/24 line ratio. 



In Figure [22] we show the [NeV] 14 /im EW vs. Sil 10 /im Strength. We see that the 
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distribution of [NeV] EW is similar in Sy l's, Sy 2's and hidden Sy l's (Table 121 13f4"f5]) . There 
is no correlation of [NeV] EW or luminosity with Sil strength (Table [TO]) . This is consistent 
with our results above that (1) the [NeV] luminosity is proportional to the radio power, and 
(2) the effect of extinction on the [NeV] 14/24 ratio is small. This suggests that either the 
torus is not very opaque at these wavelengths (14-24 /mi), or [NeV] is produced primarily 
outside the torus, e.g., in the NLR. However, we expect the high ionization lines [NeV] and 
[OIII] to be produced in similar locations around the AGN. 



Haas et al.l (120051 ) have shown that [OIII] is very orientation dependent in quasars and 
radio galaxies, suggesting that much of the [OIII] is produced in regions obscured by the 
torus. In Figure [23]we show the ratios [Offl]A5007/[OIV]26/mi and [OIII] A5007/ [NeV] 14/rni 
vs. the Sil Strength. We see a strong rank correlation for [OIII]A5007/[OIV]26/im and a 
weaker one for [OIII] A5007/ [NeV] 14/zm in the sense that the ratio is higher in objects with 
Sil emission (roughly face on) and lower in objects with Sil absorption (roughly edge on). 
We find that [OIII] is orientation dependent in our sample of Seyferts suggesting that much 
of the [OIII] emission is produced in a region which is subject to orientation dependent 
obscuration, e.g., the torus. Thus, if [NeV] is also primarily produced in regions obscured 
by the torus, then our result sugges ts that the opacity of th e torus at 14 /im is not large. 



This is consistent with the results of iBuchanan et al. I (120061 ). 



3.2.4- An Extended Hot Dust Component in Seyferts? 

High resolution mid-IR observations show extended hot dust associated with the NLR 
in NGC1068 (Galliano et al. 2005; Mason et al. 2006) and the Circinus galaxy (Packham et 
al. 2005). Figure I2"41 shows a rank correlation between the luminosities in the [NeV] 14/im 
line and the continuum determined at 6.7 /im in our Seyferts. The rank correlation is also 
present between the fluxes (Table [TO]) . In § 13.2.31 we show that the [NeV] emission is not 
strongly orientation dependent. However, in models where the 6.7 /im continuum originates 
in the hot dust in the inner torus, we would expect the 6.7 /mi continuum to be orientation 
dependent. In the context of standard uniform torus models we would expect that the 6.7 
/im emission should be high in the Sy l's, but low in the Sy 2's where it is obscured by the 
torus. Instead we see a strong rank correlation followed by both Sy l's and Sy 2's. This 
suggests that in addition to any hot dust in the torus, we are seeing an extended component 
of hot dust whose luminosity scales with the AGN luminos i ty (as measured by [NeV]). This 



might be due to dust in the NLR or a wind. iMor et al. I (120091 ) also find evidence for hot 



dust emission associated with the NLR in luminous QSOs. 
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3.3. Comparison of Seyfert Subsamples 

Here we discuss the properties of subsamples of the Seyfert population. Tabled] gives the 
median values of a range of properties for the continuum emission, Sil feature, star formation 
tracers, AGN tracers, and diagnostic line ratios. Tables [2]l9] give the results of statistical tests 
for the significance of differences in the distribution of those properties. 

The comparison of the Sy l's with the Sy 2's (both with and without a detected HBLR) 
(Table|5J) shows significant differences (p < 1%) between the Sil feature, the 20-30 /im spectral 
index, 6.7 /im continuum flux density, and the EW of the 6.2 /im PAH feature. As discussed 
in §3.2.11 the Sy l's exhibit weak Sil emission and the Sy 2's show Sil absorption. The 
other differences are consistent with being due to a higher contribution from hot dust in 
the spectra of the Sy l's which causes higher 6.7 /im continuum, bluer 20-30 /im spectral 
index, and lower EW of the 6.2 /im PAH feature. Note that although the EWs are different, 
the luminosities are similar. When we compare the Sy l's to just those Sy 2's without an 
HBLR (Table [3]) we see the same differences as between the Sy l's and all the Sy 2's, with 
the addition of a difference in the luminosity of the 6.7 /im continuum. 

The Sy 2's with and without HBLRs show differences in 20-30 /im spectral index, 6.7 
/im continuum, and EW of the 6.2 /im PAH feature and the H 2 S(l) line (Table H]). These 
differences are consistent with a higher contribution from a hot dust component in the HBLR 
Sy 2's. In addition, the Sy 2's with an HBLR are very similar to the the Sy l's; the only 
significant difference is the Sil strength (Table [5]). Thus, the Sy 2's with an HBLR are more 
similar to the Sy l's than they are to the Sy 2's without an HBLR. The differences are in the 
sense that the Sy 2's with an HBLR have a higher ratio of AGN to Starburst contribution 
to the SED. T his may contribute to the detection of the HBLR in polarized light. This is 



consistent with Alexander I (120011 ) who noted that there was a weaker contribution from the 
galaxy in the HBLR Sy 2's than the non-HBLR Sy 2's. We cannot rule out the possibility 
that the differences between the HBLR and non-HBLR Sy 2's reflect intrinsic differences in 
the AGN. 



Deo et al.l (120071 ) find that the mid-IR spectra of Sy 1.8-1.9 and are similar to those of 
Sy 2's with significant star formation. The Sy 1.8-1.9 are the smallest subgroup with only 
six objects in this study. The statistics are of course severely affected by small numbers, but 
we include the analyses for completeness. The Sy 1.8-1.9 show no significant differences with 
the Sy 2's without an HBLR (Table [6J or with an HBLR (Table W>- or the Sy l's (Table E). 
This lack of difference between the Sy 1.8-1.9 and the other Seyfert types suggests that the 
Sy 1.8-1.9 may h ave intermediate properties between the Sy l's and Sy 2's (as suggested by 



Deo et al.l (120071 )). However, it may also be due to the small number of Sy 1.8-1.9 in this 



study. 
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As noted ab ove Sy 1.8-1.9 have mid-IR spectra which are similar to those of Sy 2's 
(IDeo et all 120071 ) and HBLR Sy 2's have mid-IR spectra which are closer to those of Sy l's 
than non-HBLR Sy 2's. Thus, we also created a modified class of AGN1 consisting of the 
Sy l's plus the Sy 2's with an HBLR, and a modified class of AGN2 consisting of the Sy 2's 
without an HBLR plus the Sy 1.8-1.9. The test for differences in the properties of these two 
modified classes is shown in Table HI We compare these results (Table [9} to the comparison 
between the Sy l's and Sy 2's (Table [2]). In the modified AGN1 and AGN2 the differences in 
20-30 fim spectral index, 6.7 fim continuum flux and luminosity, and EW of the 6.2 fim PAH 
feature are stronger; the differences in Sil strength are weaker; and the modified samples 
show a new difference in the EW of the H 2 S(l) line. These findings are consistent with our 
results that the Sy 2's with an HBLR have a larger contribution from a hot dust component 
than the the Sy 2's without an HBLR. 



SUMMARY 



We present an analysis of Spitzer IRS spectroscopy of 83 active ga laxies, including 7 3 
Seyferts, 6 LINERS, and 4 HII galaxies from the exte nded 12/xm s a mple q f lRush et al. (|l993h . 
We fit the spectra using a modified version of the ISmith et all (j2007al ) PAHFIT program 
and examine selected diagnostic spectral features. 



Sy l's and Sy 2's and the lBrandl et al. I (120061 ) sample of starbursts have similar distri- 
butions of the luminosity of the 6.2/iin PAH feat ure. The EW o f the 6 .2/zm PAH feature is 
slightly larger in the Sy 2's and much larger in the lBrandl et al. I (120061 ) sample of starbursts 
than in the Sy l's due to the differences in the underlying continuum. 

We see no difference in the 6.2/11.3 PAH ratio between the different sub samples or 
any dependence of the ratio on [NeV] luminosity or Sil 10/iin strength. Thus, any effect of 
the AGN on the distribution of PAH grains must be confined to a small volume so that the 
integrated measurements studied here are not affected. 

The PAH EW is larger in sources with redder mid and far IR spectra. In addition, the 
luminosity of the 6.2/im PAH feature is rank correlated with the luminosity of the continuum 
at 60/iin. The correlation is also present in the fluxes. Both of these results suggest that 
(1) The PAH feature is a reliable tracer of star formation and (2) the there is a significant 
contribution to the heating of the cool dust by stars. The 55-90 vs. 20-30 spectral index 
plot is also a diagnostic of the relative contribution of Starburst to AGN. We see there is 
a large change in spectral index across the sample: Aa ~ 3 for both indexes. Thus, the 
contribution to the IR spectrum from the AGN and starburst components can be comparable 
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in magnitude but the relative contribution also varies widely across the sample. 

We see a rank correlation between the H 2 5(1) 17 /mi EW and two measures of the 
relative strength of the Starburst to AGN component - the 6.2/im PAH EW and the 20- 
30/im spectral index. We also see are rank correlation between the luminosities of the 6.2/im 
PAH feature and the H 2 5(1) 17 /im line. The rank correlation is also present in the fluxes. 
These results suggest that the H 2 emission is also primarily excited by star formation. 

We find rank correlations between several AGN tracers - 8 GHz vs. [NeV] 14/im [NeV] 14/im 
vs. [OIII]A5007, [GTV]26/im vs. [OIII]A5007, and [NeV] 14/im vs. 6.7/im continuum. The 
correlations are present both in the fluxes and luminosities of these AGN tracers. The rank 
correlation between [NeV] 14/tm vs. 6.7/im continuum suggests that there is a component 
to the 6.7/im continuum from hot dust which is extended, e.g., in the NLR. We find rank 
correlations of the ratios [OIII]A5007/[OIV]26/im and [OIII]A5007/[NeV] 14/im with the Sil 
strength orientation indicator. This suggests that some of the [OI H]A5007 emi s sion i n these 



Seyferts is subject to orientation dependent obscuration as found by lHaas et all (120051 ) for ra- 
dio galaxies and quasars. There is no correlation of [NeV] EW with the orientation indicator 
Sil 10/iHi strength indicating that the [NeV] emission is not strongly orientation dependent. 
This suggests with the obscu ring material (e.g . , toru s) is not very optically thick at 14/im 



consistent with the results of iBuchanan et al. I (120061 ). 



We searched for correlations between AGN and Starburst tracers (i.e., 6.2 /im PAH vs. 
[NeV] 14/im or [OIV] 26 /im or 6.7/im continuum, and [NeV] 14/im vs. 60 /im continuum). 
Although there are apparent correlations in luminosity (likely due to the dependence on 
distance), there are no significant correlations in flux, and we conclude that the AGN and 
Starburst tracers are not correlated. This is consistent with our conclusion that the relative 
strength of the AGN and Starburst components varies widely across the sample. The lack of 
a correlation between the AGN and star burst suggests that in these AGN there is no direct 
link between AGN fueling, black hole growth, and star formation. 

The density diagnostic [NeV] 14/24 /im and [SIII] 18/33 /im line ratios are consistent 
with the gas being near the low density limit, i.e., ~ 10 3 cm -3 for [NeV] and n e ~ few 
hundred cm~ 3 for [SIII]. 

The distribution of Sil 10 /tm and 18 /im strengths is consistent with the clumpy torus 



models of ISirocky et all (120081 ). 



We find a rank correlation between the [NeV] 14/im line and the continuum determined 
at 6.7 /im in our Seyferts. This suggests that in addition to any hot dust in the torus, we are 
seeing an extended component of hot dust whose luminosity scales with the AGN luminosity 
(as measured by [NeV]). This might be due to dust in the NLR or a wind. 
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We have compared the distribution of several starburst and AGN tracers between the 
different subsamples of Seyferts. Our comparison of the Sy l's with the Sy 2's without an 
HBLR shows that in addition to the strong difference in the Sil strength, the Sy 2's without 
an HBLR have 2 — 3<r differences in the 20-30 fim spectral index, 6.7 /im continuum, and the 
EW of the 6.2 /im PAH feature and [NeV] 14/xm line. The Sy 2's with an HBLR have similar 
distributions of properties to the Sy l's; the only significant difference is the Sil strength. 
Thus, the Sy 2's with an HBLR are more similar to the Sy l's than they are to the Sy 2's 
without an HBLR. The differences are in the sense that the Sy 2's with an HBLR have 
a higher ratio of AGN to Starburst contribution to the SED. Thi s may contr i bute t o the 



detection of the HBLR in polarized light. This is consistent with Alexander I (120011 ) who 
noted that there was a weaker contribution from the galaxy in the HBLR Sy 2's than the 
non-HBLR Sy 2's. 

This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. Support for this work at Bucknell University and Rochester Institute 
of Technology was provided by NASA through an award issued by JPL/Caltech. We are 
grateful to Henrik Spoon for helpful discussions and to an anonymous referee for constructive 
comments. This research made use of (1) the NASA/IPAC Extragalactic Database (NED), 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration, and (2) NASA's 
Astrophysics Data System Abstract Service. 
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Table 1. Median Values of Seyfert Properties 
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Note. - - Median values of the properties of the 4 subgroups of Seyferts. Col 1. The property. Col 2. 
Median values for the whole sample. Col 3. Median value for SI - SI. 5, and Sin. Col 4. Median value for 
SI. 8-1. 9. Col 5. Median value for Slh and Sli (hidden BLR). Col 6. Median value for S2. Luminosities 
are reported in ergs s -1 , fluxes in 10~ 14 ergs s" 1 cm -2 , and EW in nm. Medians were calculated using the 
Kaplan-Meier estimator where data include lower and upper limits (interval-censored data). Declination 
limits reduce the sample slightly for the VLA measurements. Several planned MIPS observations never 
executed, reducing slightly the number of a (55-90 jum) measurements. Non-detections of both lines also 
affect the sample size for the line ratio measurements. 
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Table 2. Comparison of Sy 1 and Sy 2 (non-HBLR + HBLR) Properties 
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Note. - The results of statistical tests to determine if the measured parameters of the 
Sy l's and Sy 2's are consistent with being drawn from the same distribution. The samples 
are defined based on optical classifications: Sy 1 includes SI - 1.5 and Sin; Sy 2 includes 
S2, Slh, Sli. Col 1. Property being tested. Col 2. Number of objects in first sample 
(Sy 1). Col 3. Number of objects in the second sample (Sy 2). Col 4. Type of statistical 
test: KS is the Kolmogorov-Smirnov test, and logrank is that test from survival analysis, 
adapted to accommodate lower and upper limits (interval-censored data). Col 5. The value 
of the test statistic; MC indicates that the probability in Col. 6 was determined by 1000 
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Monte Carlo trials. Col 6. The probability that the two samples are drawn from the same 
parent distribution of object properties; low values indicate a significant difference between 
the samples. Lower limits were removed from the sample of line ratios, but upper limits were 
included in the analysis. 
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Table 3. Comparison of Sy 1 and Sy 2 (non-HBLR) Properties 



Property 


ni 


n 2 


Test Statistic 


Value 
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Note. - The results of statistical tests to determine if the measured parameters of the 
Sy l's and Sy 2's are consistent with being drawn from the same distribution. The samples 
are defined based on optical classifications: Sy 1 includes SI— 1.5 and Sin; Sy 2 includes S2 
only. Col I. Property being tested. Col 2. Number of objects in first sample (Sy 1). Col 
3. Number of objects in the second sample (Sy 2). Col 4. Type of statistical test: KS is 
the Kolmogorov-Smirnov test, and logrank is that test from survival analysis, adapted to 
accommodate lower and aupper limits (interval-censored data). Col 5. The value of the 
test statistic. MC indicates that the probability in Col. 6 was determined by 1000 Monte 
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Carlo trials. Col 6. The probability that the two samples are drawn from the same parent 
distribution of object properties; low values indicate a significant difference between the 
samples. Lower limits were removed from the sample of line ratios, but upper limits were 
included in the analysis. 
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Table 4. Comparison of HBLR and non-HBLR Sy 2 Properties 



Property 
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Note. — The results of statistical tests to determine if the measured parameters of hidden 
Sy l's and standard Sy 2's are consistent with being drawn from the same distribution. The 
samples are defined based on optical classifications, with HBLR comprising Sli and Slh, and 
non-HBLR comprising only S2. Col 1. Property being tested. Col 2. Number of objects in 
first sample (Sy 1). Col 3. Number of objects in the second sample (Sy 2). Col 4. Type of 
statistical test: KS is the Kolmogorov-Smirnov test, and logrank is that test from survival 
analysis, adapted to accommodate lower and upper limits (interval-censored data). Col 5. 
The value of the test statistic. MC indicates that the probability in Col. 6 was determined 
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by 1000 Monte Carlo trials. Col 6. The probability that the two samples are drawn from 
the same parent distribution of object properties; low values indicate a significant difference 
between the samples. Lower limits were removed from the sample of line ratios, but upper 
limits were included in the analysis. 
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Table 5. Comparison of Sy 1 and HBLR Sy 2 Properties 
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0.19 


76.6 


SJQ.7 urn) 


27 


17 


KS 


0.22 


61.7 


log z/L^(6.7 /im) 


27 


17 


KS 


0.31 


20.9 


IRAS log z/L„(60 urn) 


27 


17 


KS 


0.22 


67.0 


Sil 10 /im 


27 


17 


KS 


0.51 


0.8 


log uLJ8A GHz) 


24 


16 


logrank 


MC 


24.8 


EW(PAH 6.2 urn) 


27 


17 


logrank 


MC 


70.8 


log L(PAH 6.2 /im) 


27 


17 


logrank 


MC 


95.4 


EW(H 2 S(l) A17/im) 


27 


17 


logrank 


MC 


32.6 


log L(H 2 S(l) A17/im) 


27 


17 


logrank 


MC 


61.2 


log L(S III A18 p) 


27 


17 


logrank 


MC 


60.6 


EW(Ne V A14 /im) 


27 


17 


logrank 


MC 


20.0 


log L(Ne V A14 /im) 


27 


17 


logrank 


MC 


18.4 


log L(0 IV A26 /im) 


27 


17 


logrank 


MC 


26.8 


Ne V 14 /im / 24 /im 


20 


13 


logrank 


MC 


72.8 


S III 18 /im / 33 /im 


19 


12 


logrank 


MC 


66.8 


PAH 6.2 /im / 11.3 /im 


26 


16 


logrank 


MC 


90.2 



Note. — The results of statistical tests to determine if the measured parameters of Sy l's 
and Sy 2's with HBLR are consistent with being drawn from the same distribution. The 
samples are defined based on optical classifications. Col 1. Property being tested. Col 2. 
Number of objects in first sample (Sy 1). Col 3. Number of objects in the second sample 
(Sy 2). Col 4. Type of statistical test: KS is the Kolmogorov-Smirnov test, and logrank is 
that test from survival analysis, adapted to accommodate lower and upper limits (interval- 
censored data). Col 5. The value of the test statistic. MC indicates that the probability 
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in Col. 6 was determined by 1000 Monte Carlo trials. Col 6. The probability that the two 
samples are drawn from the same parent distribution of object properties; low values indicate 
a significant difference between the samples. Lower limits were removed from the sample of 
line ratios, but upper limits were included in the analysis. 
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Table 6. Comparison of Sy 1.8-1.9 and non-HBLR Sy 2 Properties 



Property 


ni 


n 2 


Test Statistic 


Value 


Two-sided significance (%) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




Distance 


6 


20 


KS 


0.35 




52.8 


a(20 - 30 fan) 


6 


20 


KS 


0.38 




40.9 


a(55 — 90 jum) 


6 


19 


KS 


0.26 




82.2 


5^(6.7 /im) 


6 


20 


KS 


0.33 




59.3 


log uL u (6.7 /im) 


6 


20 


KS 


0.35 




52.8 


IRAS log uL u (Q0 /im) 


6 


20 


KS 


0.20 




97.2 


Sil 10 /xm 


6 


20 


KS 


0.33 




59.3 


log i/L„(8.4 GHz) 


6 


19 


logrank 


MC 




68.6 


EW(PAH 6.2 /im) 


6 


20 


logrank 


MC 




59.4 


log L(PAH 6.2 /im) 


6 


20 


logrank 


MC 




91.2 


EW(H 2 S(l) A17/im) 


6 


20 


logrank 


MC 




31.6 


log L(H 2 S(l) A17/im) 


6 


20 


logrank 


MC 




25.2 


log L(S III A18 /im) 


6 


20 


logrank 


MC 




32.2 


EW(Ne V A14 /im) 


6 


20 


logrank 


MC 




38.2 


log L(Ne V A14 /im) 


6 


20 


logrank 


MC 




49.2 


log L(0 IV A26 /im) 


6 


20 


logrank 


MC 




72.6 


Ne V 14 /im / 24 /im 


4 


13 


KS 


0.29 




91.9 


S III 18 /im / 33 /im 


5 


16 


logrank 


MC 




13.2 


PAH 6.2 /im / 11.3 /im 


6 


20 


logrank 


MC 




37.6 



Note. — The results of statistical tests to determine if the measured parameters of Sy 1.8- 
1.9 and Sy 2's without HBLR are consistent with being drawn from the same distribution. 
The samples are defined based on optical classifications. Col 1. Property being tested. Col 
2. Number of objects in first sample (Sy 1). Col 3. Number of objects in the second sample 
(Sy 2). Col 4. Type of statistical test: KS is the Kolmogorov-Smirnov test, and logrank is 
that test from survival analysis, adapted to accommodate lower and upper limits (interval- 
censored data). Col 5. The value of the test statistic. MC indicates that the probability 
in Col. 6 was determined by 1000 Monte Carlo trials. Col 6. The probability that the two 
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samples are drawn from the same parent distribution of object properties; low values indicate 
a significant difference between the samples. Lower limits were removed from the sample of 
line ratios, but upper limits were included in the analysis. 
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Table 7. Comparison of Sy 1.8-1.9 and HBLR Sy 2 Properties 



Property 


ni 


n 2 


Test Statistic 


Value 


Two-sided significance (%) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




Distance 


6 


17 


KS 


0.26 




84.7 


a(20 - 30 fan) 


6 


17 


KS 


0.44 




25.0 


a(55 — 90 jum) 


6 


17 


KS 


0.66 




2.4 


5^(6.7 /im) 


6 


17 


KS 


0.44 




25.0 


log uL u (6.7 /im) 


6 


17 


KS 


0.27 




79.3 


IRAS log uL u (Q0 /im) 


6 


17 


KS 


0.44 




25.0 


Sil 10 /xm 


6 


17 


KS 


0.35 




54.1 


log i/L„(8.4 GHz) 


6 


16 


logrank 


MC 




17.6 


EW(PAH 6.2 /im) 


6 


17 


logrank 


MC 




3.2 


log L(PAH 6.2 /im) 


6 


17 


logrank 


MC 




54.4 


EW(H 2 S(l) A17/im) 


6 


17 


logrank 


MC 




7.4 


log L(H 2 S(l) A17/im) 


6 


17 


logrank 


MC 




81.2 


log L(S III A18 /im) 


6 


17 


logrank 


MC 




67.4 


EW(Ne V A14 /im) 


6 


17 


logrank 


MC 




72.0 


log L(Ne V A14 /im) 


6 


17 


logrank 


MC 




38.2 


log L(0 IV A26 /im) 


6 


17 


logrank 


MC 




44.8 


Ne V 14 /im / 24 /im 


4 


13 


logrank 


MC 




38.8 


S III 18 /im / 33 /im 


5 


12 


logrank 


MC 




39.0 


PAH 6.2 /im / 11.3 /im 


6 


16 


logrank 


MC 




3.4 



Note. - The results of statistical tests to determine if the measured parameters of Sy 
1.8-1.9 and Sy 2's with HBLR are consistent with being drawn from the same distribution. 
The samples are defined based on optical classifications. Col 1. Property being tested. Col 
2. Number of objects in first sample (Sy 1). Col 3. Number of objects in the second sample 
(Sy 2). Col 4. Type of statistical test: KS is the Kolmogorov-Smirnov test, and logrank is 
that test from survival analysis, adapted to accommodate lower and upper limits (interval- 
censored data). Col 5. The value of the test statistic. MC indicates that the probability 
in Col. 6 was determined by 1000 Monte Carlo trials. Col 6. The probability that the two 



-40- 



samples are drawn from the same parent distribution of object properties; low values indicate 
a significant difference between the samples. Lower limits were removed from the sample of 
line ratios, but upper limits were included in the analysis. 
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Table 8. Comparison of Sy 1.8-1.9 and Sy 1 Properties 



Property 


ni 


n 2 


Test Statistic 


Value 


Two-sided significance (%) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




Distance 


6 


27 


KS 


0.17 




99.6 


a(20 - 30 fan) 


6 


27 


KS 


0.65 




1.8 


a(55 — 90 jum) 


6 


26 


KS 


0.54 




7.9 


5^(6.7 /im) 


6 


27 


KS 


0.46 




18.5 


log uL u (6.7 /im) 


6 


27 


KS 


0.31 




62.2 


IRAS log uL u (Q0 /im) 


6 


27 


KS 


0.39 




44.8 


Sil 10 /xm 


6 


27 


KS 


0.43 




26.6 


log i/L„(8.4 GHz) 


6 


24 


logrank 


MC 




16.8 


EW(PAH 6.2 /im) 


6 


27 


logrank 


MC 




1.2 


log L(PAH 6.2 /im) 


6 


27 


logrank 


MC 




59.8 


EW(H 2 S(l) A17/im) 


6 


27 


logrank 


MC 




30.4 


log L(H 2 S(l) A17/im) 


6 


27 


logrank 


MC 




57.0 


log L(S III A18 /im) 


6 


27 


logrank 


MC 




95.0 


EW(Ne V A14 /im) 


6 


27 


logrank 


MC 




39.6 


log L(Ne V A14 /im) 


6 


27 


logrank 


MC 




84.4 


log L(0 IV A26 /im) 


6 


27 


logrank 


MC 




65.6 


Ne V 14 /im / 24 /im 


4 


20 


logrank 


MC 




48.8 


S III 18 /im / 33 /im 


5 


19 


logrank 


MC 




35.2 


PAH 6.2 /im / 11.3 /im 


6 


26 


logrank 


MC 




20.4 



Note. - The results of statistical tests to determine if the measured parameters of Sy 
1.8-1.9 and Sy l's are consistent with being drawn from the same distribution. The samples 
are defined based on optical classifications. Col I. Property being tested. Col 2. Number 
of objects in first sample (Sy 1). Col 3. Number of objects in the second sample (Sy 2). 
Col 4. Type of statistical test: KS is the Kolmogorov-Smirnov test, and logrank is that test 
from survival analysis, adapted to accommodate lower and upper limits (interval-censored 
data). Col 5. The value of the test statistic. MC indicates that the probability in Col. 6 
was determined by 1000 Monte Carlo trials. Col 6. The probability that the two samples are 
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drawn from the same parent distribution of object properties; low values indicate a significant 
difference between the samples. Lower limits were removed from the sample of line ratios, 
but upper limits were included in the analysis. 
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Table 9. Comparison of (Sy 1 + HBLR Sy 2) and (Sy 1.8-1.9 + non-HBLR Sy 2) 

Properties 



Property 


ni 


™2 


Test Statistic 


Value 


Two-sided significance (%) 


I 1 ) 


( 2 ) 


{Q\ 
( 3 ) 


\V 


\P) 


(6) 


Distance 


44 


26 


KS 


0.18 


65.2 


a(20 - 30 jum) 


44 


26 


KS 


0.55 


< 0.1 


a(55 — 90 /im) 


43 


25 


KS 


0.39 


1.1 


S u (6.7 /im) 


44 


26 


KS 


0.55 


< 0.1 


log pL u {Q.7 /im) 


44 


26 


KS 


0.44 


0.2 


IRAS log uL u (Q0 /an) 


44 


26 


KS 


0.26 


23.7 


Sil 10 /im 


44 


26 


KS 


0.35 


3.5 


log uL u (8A GHz) 


40 


25 


logrank 


MC 


5.6 


EW(PAH 6.2 /im) 


44 


26 


logrank 


MC 


0.2 


log L(PAH 6.2 /im) 


44 


26 


logrank 


MC 


28.2 


EW(H 2 S(l) A17/im) 


44 


26 


logrank 


MC 


0.6 


log L(H 2 S(l) A17/im) 


44 


26 


logrank 


MC 


19.6 


log L(S III A18 /xm) 


44 


26 


logrank 


MC 


14.4 


EW(Ne V A14 /im) 


44 


26 


logrank 


MC 


28.4 


log L(Ne V A14 /im) 


44 


26 


logrank 


MC 


13.2 


log L(0 IV A26 /im) 


44 


26 


logrank 


MC 


34.0 


Ne V 14 /im / 24 /im 


33 


17 


logrank 


MC 


8.2 


S III 18 /im / 33 /im 


31 


21 


logrank 


MC 


42.0 


PAH 6.2 /im / 11.3 /im 


42 


26 


logrank 


MC 


3.8 



Note. - The results of statistical tests to determine if the measured parameters of the 
combined groups (Sy 1 + HBLR Sy 2) and (Sy 1.8-1.9 and non-HBLR Sy 2) are consistent 
with being drawn from the same distribution. The samples are defined based on optical 
classifications. Col 1. Property being tested. Col 2. Number of objects in first sample 
(Sy 1). Col 3. Number of objects in the second sample (Sy 2). Col 4. Type of statistical 
test: KS is the Kolmogorov-Smirnov test, and logrank is that test from survival analysis, 
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adapted to accommodate lower and upper limits (interval-censored data). Col 5. The value 
of the test statistic. MC indicates that the probability in Col. 6 was determined by 1000 
Monte Carlo trials. Col 6. The probability that the two samples are drawn from the same 
parent distribution of object properties; low values indicate a significant difference between 
the samples. Lower limits were removed from the sample of line ratios, but upper limits were 
included in the analysis. 
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Table 10. Tests for Correlations in Seyfert Properties 



First Pronertv 


Serond Pronertv 


z 


v(%) 


(1) 


(2) 


(3) 


(4) 


EW(PAH 6.2 /im) 


Distance 


-2.058 


4.0 


EW(PAH 6.2 /im) 


Sil 10 /im 


-2.311 


2.1 


EW(PAH 6.2 /im) 


a(20 - 30 /im) 


7.691 


< 0.1 


EW(PAH 6.2 //m) 


a (55 — 90 /xm) 


5.663 


< 0.1 


EW(PAH 6.2 /^m) 


EW(H 2 S(l) A17/im) 


4.888 


< 0.1 


EW(PAH 6.2 /im) 


EW(Ne V A14 /im) 


2.220 


2.6 


EW(PAH 6.2 /im) 


log L(Ne V A 14 /im) 


-0.668 


50.4 


PAH 6.2 /im Flux 


Ne V A14 /im Flux 


1.933 


5.3 


log L(PAH 6.2 /im) 


log L(Ne V A14 /im) 


4.205 


< 0.1 


PAH 6.2 /im Flux 


H 2 S(l) A17/im Flux 


4.969 


< 0.1 


log L(PAH 6.2 /im) 


log L(H 2 S(l) A17/im) 


5.333 


< 0.1 


log L(PAH 6.2 /im) 


a(20 - 30 /im) 


2.671 


0.8 


log L(PAH 6.2 /im) 


a(55 — 90 /im) 


1.370 


17.1 


PAH 6.2 /im Flux 


IRAS ^(60 /im) 


7.260 


< 0.1 


log L(PAH 6.2 /im) 


IRAS log 1/2^(60 /im) 


7.491 


< 0.1 


PAH 6.2 /im Flux 


IV A26 /im Flux 


1.029 


30.3 


log L(PAH 6.2 /im) 


log L(0 IV A26 /im) 


2.744 


0.6 


a (20 - 30 /im) 


a (55 — 90 /im) 


4.288 


< 0.1 


a (20 - 30 /im) 


EW(H 2 S(l) A17/im) 


4.292 


< 0.1 


a (20 - 30 /im) 


IRAS log vL„{60 /im) 


2.788 


0.5 


a (20 - 30 /im) 


EW(Ne V A14 /im) 


2.933 


0.3 


a (20 - 30 /im) 


log L(Ne V A 14 /im) 


-0.317 


75.2 


a (20 - 30 /im) 


log L(S III A18 /im) 


2.861 


0.4 


a (20 - 30 /im) 


Sil 10 /im 


-4.674 


< 0.1 


S v (6.7 /im) 


Ne V A14 /im Flux 


3.228 


0.1 


log z/L^(6.7 /im) 


log L(Ne V A14 /im) 


5.136 


< 0.1 


^(6.7 /im) 


PAH 6.2 /im Flux 


-0.690 


49.0 


log z/Lj,(6.7 /im) 


log L(PAH 6.2 /im) 


2.884 


0.4 


EW(Ne V A14 /im) 


Sil 10 /im 


-2.997 


0.3 


Ne V A14 /im Flux 


^(8.4 GHz) 


3.118 


0.2 



Table 10 — Continued 



First Property 


Second Property 


z 


p(%) 


(1) 


(2) 


(3) 


(4) 


log L(Ne V A14 /im) 


log z/L„(8.4 GHz) 


4.371 


< 0.1 


Ne V A14 /im Flux 


IRAS S„(60 /im) 


2.724 


0.6 


log L(Ne V A14 /im) 


IRAS log z/L^(60 /im) 


4.808 


< 0.1 


log L(Ne V A14 /im) 


Sil 10 /im 


-0.410 


68.1 


Ne V 14 /im / 24 /im 


S III 18 /im / 33 /im 


-0.171 


86.4 


Ne V 14 /im / 24 /im 


Sil 10 /tm 


-0.914 


36.1 


PAH 6.2 /im / 11.3 /im 


log L(Ne V A14 /im) 


0.501 


61.7 


PAH 6.2 /im / 11.3 /im 


Sil 10 /im 


-0.006 


99.5 


X-ray Flux (2-10 keV) 


Ne V A 14 /im Flux 


1.685 


9.2 


III Flux 


Ne V A14 /im Flux 


1.811 


7.0 


Ne V 14 /im / [Nell] 


EW(PAH 6.2 /im) 


-5.261 


< 0.1 


III / Ne V (14 /im) 


Sil 10 /im 


4.197 


< 0.1 


III / IV (26 /im) 


Sil 10 /im 


4.121 


< 0.1 


log L(0 IV A26 /im) 


log L(0 III ) 


4.847 


< 0.1 


IV A26 /xm Flux 


III Flux 


2.564 


1.0 


log L(Ne V A 14 /im) 


log L(0 III ) 


4.655 


< 0.1 


Ne V A 14 /im Flux 


III Flux 


1.811 


7.0 


log L(Ne V A 14 /im) 


Sil 10 /tm 


-0.410 


68.1 


A" H (cm" 2 ) 


Sil 10 /im 


-2.372 


1.8 



Note. — The results of non-parametric rank tests to look for cor- 
relations among the properties of the combined sample. Col 1 & 
2. Properties tested for correlation. Col 3. Kendall's generalized 
rank correlation coefficient, adjusted to account for upper and lower 
limits. Reported is the z-s c ore, f ollowing the method outlined by 



Isobe. Feigelson. fc Nelson I (119861 ). Col 4. Two-sided significance, 



giving the probability of no correlation; small values (p < 1%) indi- 
cate evidence for correlation. 
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Brandl et al. (2006) Starbursts 



0.5 



1.0 



1.5 2.0 2.5 3.0 
log EW(PAH 6.2 ^m) (nm) 



3.5 



4.0 



4.5 



Fig. 1. — The distribution of 6.2 /mi PAH EW. From the top, the panels are (1) Sy 1 
and Sy 2; (2) Sy 1 (Sl.0-1.9, Sin), (3) Sy 2 (S2, Slh, & Sli), and starbursts from the 
Brandl et al. I ((20061) sample. Note that the Brandl et al. sample shares two objects in 
common with our sample, NGC 1097 and NGC1365. The starburst measurements are derived 
from a PAHFIT spectral decomposition of the Brandl et al. data. 
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Fig. 2. — The distribution of 6.2 /zm PAH luminosity. From the top, the panels are (1) 
Sy 1 and Sy 2; (2 ) Sy 1 (SI. 0-1. 9, Sin), (3) Sy 2 (SI, Slh, & Sli), and starbursts from 



the iBrandl et al. I (120061 ) sample. Note that the Brandl et al. sample shares two objects 
in common with our sample, NGC 1097 and NGC1365. The starburst measurements are 
derived from a PAHFIT spectral decomposition of the Brandl et al. data. 
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Fig. 3. — The ratio of 6.2/11.3 /xm PAH strength as a function of [NeV] luminosity. Open 
circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, 
partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are LINERs, 
and filled stars are Starburst/HII galaxies. Limits are indicated by appropriately directed 
arrows. The single lower limit in the PAH ratio is MRK 704, for which the 11.3 /zm feature 
is not detected. 
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Fig. 4. — Ratio of the 6.2 and 11.3 fj,m PAH lines vs. Sil 10 fj,m strength. For the Sil 10 /zm 
strength, positive means emission. (Left) the data. (Right) a blow up of the region around a 
Sil strength of zero. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, 
filled circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares 
labeled L are LINERs, and filled stars are Starburst/HII galaxies. Limits are indicated by 
appropriately directed arrows. 
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Fig. 5. — 6.2 fim PAH luminosity vs. IRAS uL u (60 //m). Open circles are SI, SI. 5, and Sin, 
gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden 
Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII 
galaxies. 
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-4 -3 -2 -1 

a(55 - 90 /im) 

Fig. 6.— 6.2 PAH EW vs. slope of the MIPS SED. The slope of the MIPS SED indicates 
the relative importance of cold vs. warm dust to the FIR SED. Positive slope is a red MFIR 
SED. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles 
are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are 
LINERs, and filled stars are Starburst/HII galaxies. 
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-3-2-10 1 2 3 

a(20 - 30 yuim) 

Fig. 7. — 6.2 fim PAH EW vs. spectral index of the continuum between 20 and 30 /xm. 
Positive slope is a red MFIR SED. Open circles are SI, SI. 5, and Sin, gray-filled circles are 
SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and 
Sli), squares labeled L are LINERs, and filled stars are Starburst/HIl galaxies. 
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Fig. 8. — Luminosity of the PAH 6.2 /im feature vs. spectral index of the continuum between 
20 and 30 /im. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled 
circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled 
L are LINERs, and filled stars are Starburst/HII galaxies. 
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Fig. 9. — MIPS SED slope vs. IRS 20-30 /im spectral index. Sources with blue or flat MIPS 
SEDs tend to have flat 20-30 SEDs. Sources with red MIPS SEDs tend to have red 20-30 
SEDs with a broad distribution. Open circles are SI, SI. 5, and Sin, gray-filled circles are 
SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and 
Sli), squares labeled L are LINERs, and filled stars are Starburst/HII galaxies. 
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EW(H 2 S(l) A17^m) (nm) 

Fig. 10.— 6.2 /im PAH EW vs. H 2 S(l) 17 fjm EW. Open circles are SI, S1.5, and Sin, 
gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden 
Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII 
galaxies. Limits are indicated by appropriately directed arrows. 
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Fig. 11. — IRS 20-30 //m spectral index vs. H 2 S(l) 17/xm EW. We see a trend for sources 
with redder 0:20-30 to have larger molecular line EW. This is consistent with molecular 
hydrogen emission and a red SED both being signatures of star formation. Open circles are 
SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled 
circles are "hidden Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars 
are Starburst/HII galaxies. Limits are indicated by appropriately directed arrows. 
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Fig. 12. — Luminosity of the [NeV] 14/xm line vs. the luminosity of the 8 GHz radio contin- 
uum. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles 
are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are 
LINERs, and filled stars are Starburst/HII galaxies. Limits are indicated by appropriately 
directed arrows. 
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Fig. 13. — (Left). Luminosity of the [NeV] 14/zm line vs. the luminosity of the [OIII]A5007 
line. (Right). Luminosity of the [OIV] 26/im line vs. the luminosity of the [OIII]A5007 line. 
Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 
2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are LINERs, 
and filled stars are Starburst/HII galaxies. Limits are indicated by appropriately directed 
arrows. 
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Fig. 14. — 6.2 fim PAH luminosity vs. [NeV] 14 /iin luminosity. Open circles are SI, SI. 5, 
and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles 
are "hidden Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are 
Starburst/HII galaxies. Limits are indicated by appropriately directed arrows. 
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Fig. 15. — 6.2 yum PAH feature flux vs. [NeV] 14 /itn line flux. This is an AGN vs Starburst 
relation. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles 
are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are 
LINERs, and filled stars are Starburst/HII galaxies. Limits are indicated by appropriately 
directed arrows. 
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Fig. 16.— The ratio [NeV] 14/zm/ [Nell] 12/xm vs. 6.2//m PAH EW. Open circles are SI, SI. 5, 
and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles 
are "hidden Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are 
Starburst/HII galaxies. 
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Fig. 17. — 6.2 [im PAH EQW vs. Sil 10 /iin strength. For the Sil 10 strength, positive 
means emission. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, 
filled circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares 
labeled L are LINERs, and filled stars are Starburst/HII galaxies. Limits are indicated by 
appropriately directed arrows. 
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Fig. 18. — Sil 18 um Strength vs. Sil 10 /im strength. For the Sil strength, positive means 
emission. (Left) All the points. Blow up of inset region is on the right. (Right) Blow up of 
region near Sil Strength of zero. Open circles are SI, SI. 5, and Sin, gray-filled circles are 
SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden Seyferts" (Slh and 
Sli), squares labeled L are LINERs, and filled stars are Starburst/HII galaxies. Limits are 
indicated by appropriately directed arrows. 



-65 - 



0.6 



0.4 



0.2 



to 0.0 



-0.2 



-0.4 



-0.6 



q = 



x= 100 




□ 



'''"A. 



;A N = 2 



S. 1 



"A N = 3 



*V. n = 5 



A—A N = 4 

A 



■• A. A 



A' 



A 



N= 10 



'.•"•■A' N = 6 
N = 8 



-2.0 



-1.5 



-1.0 



-0.5 



0.0 



SillO 



0.5 



1.0 



Fig. 19. — Sil 18 /im Strength vs. Sil 10 /im strength. For the S il strength, positive 
means emission. Theoretical curves for spherical obscuration models (jSirocky et al.l 120081 ) 
are shown. Dotted lines are clumpy models and dashed lines are smooth models. Model 
parameters are annotated: q is the power law index of the radial density distribution (here, 
flat); Y is the ratio of the outer radius to the inner radius of the obscuring region (spherical 
models, the inner radius is set by dust sublimation); r refers to r v , the total opacity through 
the spherical shell (spherical models); and N refers to the average number of clumps along 
the sight-line (clumpy models). Spherical models range over Y = 100, 200, 300, 400, 500, 600, 
800, & 1000; Y = 30 for the clumpy models. The clumpy models are further distinguished 
by the (unlabeled) opacity through individual clumps: r v = 10, 20, 30, 40, 60, & 80, with 
lower opacities oriented to the upper right on this diagram. 
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Fig. 20. — [NeV] 14//m/24//m vs. [SIII] 18//m/33/im. Open circles are SI, SI. 5, and Sin, 
gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden 
Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII 
galaxies. Limits are indicated by appropriately directed arrows. 
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Fig. 21. — [NeV] 14/im/24/zm vs. Sil 10 /im strength. Open circles are SI, SI. 5, and Sin, 
gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden 
Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII 
galaxies. Limits are indicated by appropriately directed arrows. 




Fig. 22. — [NeV] 14/xm EW vs. Sil 10 /im strength. Open circles are SI, SI. 5, and Sin, 
gray-filled circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden 
Seyferts" (Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII 
galaxies. Limits are indicated by appropriately directed arrows. 
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Fig. 23.— (Left) The ratio [OIII]5007/[OIV]26//m vs. Sil 10 /im strength. (Right) The ratio 
[OIII]5007/[NeV]14/zm vs. Sil 10 /im strength. Open circles are SI, SI. 5, and Sin, gray-filled 
circles are SI. 8 and SI. 9, filled circles are Sy 2, partially filled circles are "hidden Seyferts" 
(Slh and Sli), squares labeled L are LINERs, and filled stars are Starburst/HII galaxies. 
Limits are indicated by appropriately directed arrows. 
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Fig. 24. — Luminosity of the [NeV] 14//m line vs. the luminosity of the continuum at 6.7 
/im. Open circles are SI, SI. 5, and Sin, gray-filled circles are SI. 8 and SI. 9, filled circles 
are Sy 2, partially filled circles are "hidden Seyferts" (Slh and Sli), squares labeled L are 
LINERs, and filled stars are Starburst/HII galaxies. Limits are indicated by appropriately 
directed arrows. 



